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Abstract 

; The baryon differential spectrum of the baryon decay A;, — t- h^civi will 



O 



be measured in detail at LHCb. We obtain new results on the form factors 
in the heavy quark expansion of Heavy Quark Effective Theory that can be 



> 

o 

■ useful in the interpretation of the data. We formulate a sum rule for the 



elastic subleading form factor A{w) at order l/mq, that originates from the 
Lagrangian perturbation Ckin- In the sum rule appear only the intermediate 



states {j^ , J^) = ^0''",^''^^, entering also in the l/niQ correction to the 
axial form factor Gi{w), that contributes to the differential rate at zero recoil 
^ ■ w = 1. This result, together with another sum rule in the forward direction 

I for |Gi(l)p, allows us to obtain a lower bound for the correction at zero recoil 

—Sy^2 in terms of the derivative ^i(l) and the slope p\ and curvature cr^ of 
the elastic Isgur-Wise function ^\{w). Another theoretical implication is that 
^'(1) must vanish for some relation between p\ and a\, as well as for p\ — t- 0, 
establishing a non-trivial correlation between the leading IW function ^a(^^) 
and the subleading one A{w). A phenomenological estimation of these two 
functions allows to obtain a lower bound on —6) / , . 
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1 Introduction 



The present paper concerns the apphcation, within Heavy Quark Effective The- 
ory (HQET), of the formahsm used to study the decay B — t- D^*^ii'e for the 
subleading corrections in l/mg due to the Lagrangian perturbation in the baryon 
transition A;, — )■ Aciue- In the present paper, as in [1], one important ingredient is 
the consideration of the non-forward amphtude Ah{vi) — Ac{v') — ?■ Ah{vf), allowing 
for general four- velocities Vi, Vf, v' , first proposed by Uraltsev in the meson case at 
leading order of the heavy mass expansion [2]. 

In the meson case, at leading order, the Bjorken Sum Rule (SR) [3][1] gives the 
lower bound for the derivative of the meson elastic IW function at zero recoil w = 1, 

= —.^'(1) > J- In [2], Uraltsev used the non-forward amplitude to derive a new 
SR that, combined with Bjorken's, gave the much stronger lower bound > |. In 
[5], we did develop a manifestly covariant formalism within the Operator Product 
Expansion (OPE) and the non-forward amplitude, using the whole tower of heavy 
meson states [Ej. We did recover Bjorken and Uraltsev SR plus a general class 
of SR that also allow us to bound higher derivatives of the IW function [7]. In 
particular, we found a bound on the curvature in terms of the slope p^, namely 

Recently, at leading order, we have extended the SR method to the elastic leading 
IW function ^aI^) of the baryonic transition Af, — )■ A^fz/^ [9]. We have recovered 
the lower bound for the slope p\ = —,^^(1) > [10] and generalized it, as we briefly 
summarize in Section 3. 

The shape of the differential spectrum of the decay Af, — > Acii'e depends both on 
these results at leading order, as well as on the ©(l/mg) corrections, which are the 
main subject of the present paper. The mode Af, — )■ Acii^e has a large branching ratio 
already measured at the Tevatron, of about 5 %, a large fraction of the inclusive 
semileptonic decay BR{Ab — > Acii'i + anything) = 10%. At LHCb, precise data on 
the decay A;, — )■ A^fz/^, in particular on the differential rate, will be obtained in near 
future. 

The paper is organized as follows. Section 2 makes explicit the notations on the 
form factors and the heavy quark expansion. In Section 3, we summarize results 
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on the leading IW function $,Aiw) for baryons. Section 4 gives a short derivation 
of the relevant SR for the Lagrangian elastic I/ttlq axial form factor A{w) in the 
transition A;, — A^, by starting from the definition of the subleading Lagrangian 
form factors. In Section 5, we deduce another SR for the 0{l/mQ) correction at 
zero recoil of the axial form factor Gi(w), denoted — 5, , 2 , that contributes to the 

^ ' ' 1 / THq ' 

differential rate a.t w = 1. We show that the same states with light cloud = O"*" 
contributing to this latter SR, contribute also to another SR for A{w). In Section 
6, using Schwarz inequality, we obtain a lower bound on the correction —6\'l^\ in 



terms of the subleading elastic form factor A{w) and the leading elastic IW function 
^a(w^)- In Section 7, we establish a correlation between A'{1) and the shape of 
the leading IW function and in Section 8 we summarize the physical parameters 
involved in the bound. In Section 9, following Leibovich and Stewart [11], we bound 
one of these parameters. In Section 10, we summarize theoretical results on leading 
and subleading IW functions that follow from HQET and from QCD Sum Rules 
(QCDSR). In Section 11, we perform a calculation of A'{1) in the quark model. 
Finally, in Section 12, we give our numerical results on the lower bound on —5, , 2 , 



and in Section 13, we conclude. 

2 Baryon leading and subleading IW functions 

We start defining the notations for the matrix elements and the form factors follow- 
ing Falk and Neubert [12] : 



{A,{v\ s') I cYb I A,(^;, s)) = ua^v, s') [^17^^ + + F^v'^] ua^v, s) , (1) 
{A,{v', s') I cYl5b I Ab{v, s)) = UA^v', s') [G.^ + G2V>' + G^v'^^] j^UA.iv, s) (2) 



where ua(v, s) are spinors normalized by ua{v, s)ua{v, s') = 2mA5s,s'. 

In the heavy quark mass expansion, up to the order l/rriQ, these form factors 
are given in terms of two functions Bi{w) and B2{w) of w = v ■ v' [13] : 



■Q 



■Q 




(3) 



(4) 
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F^{w) = G2{w) = ^B2{w) , (5) 



2m, 



c 



F,{w) = -G^{w) = ^B2{w) (6) 



where 



w — 1 - 

i?iH = ^- AaH + AH, (7) 

w + 1 

B2{w) = ^ A ■ (8) 

w + 1 

For a generic heavy quark current J = h^iVh^ (F is any Dirac matrix), the elastic 
leading IW function is defined by: 

{k,{v\ s') I h^)Vh^^ I K,{v, s)) = a(w^)^A.(t/, s')Vu!,^{v, s) . (9) 

The terms proportional to A^a('U^) come from the l/mq Current-type pertur- 
bations. Note that the situation for mesons is different. In this case, there are 
two types of current perturbations, namely, one proportional to AC,{w) and another 
independent function C,3{uj) [H]. 

The function A{w) comes from the kinetic part Ckin,v of the Lagrangian pertur- 
bation: 

^eff = ^0 + ^kin,v + ^mag,v + 0{l/mQ) , (10) 
A.n,. = ^Oi2,. , O^Sl = h^^H^D^W\ (11) 

'-mag,!! 2mQ ^^"'^'^ ' rnag,v 2 v 1^'^ v ■ l"^"^/ 

For a current J = hy'Th^, one has: 

{A,{v', s') I z f d'xT{J{0), Oitix)} I Mv, s)) = A{w)uaJua, . (13) 



Notice that in the transition A^ — )■ K^ii'i under consideration, the magnetic part 
^mag,v of the l/mq pcrturbatiou to the Lagrangian does not contribute [13]. This 
is unlike the case of the meson ground state = | which involves instead three 
independent Lagrangian perturbation functions, coming from both Ckm,v and Cmag,v 
and denoted by Xj(^) (i = 1)2,3) [1] 



An important physical remark must be emphasized here, namely, the simplicity 
of the baryon transition A;, — )■ h^c^^vi concerning the l/mg corrections. Indeed, 
we have a single Current-tjpe perturbation proportional to A^\{w) and a single 
Lagrangian-tjpe form factor A{w) originating from the kinetic part Ckin,v This 
simplification compared to the meson case comes merely from the fact that the light 
cloud in the baryon ground state Aq has = 0~^. This is in sharp contrast with the 
the meson ground state decay B D{D*)ih'i, where the light cloud has = ^ ■ 

The subleading baryon form factor A{w) coming from Ckm.v satisfies, due to 
vector current conservation, the condition at zero recoil |12j : 

A{1) = . (14) 

The differential rate for A;, — Aciue in the neighborhood of the zero recoil point 
w = 1 depends only on the form factor Gi{w). One gets, in the limit w —> 1: 

vJ-i '"'^'"'7^'""'' " ^ij^^'"^- - '"A.)'iG.(m' (15) 

where Gi(l) has only corrections to order I/itiq as made explicit in formula (5.2) 
of [12]: 

We will come back to the detail of these corrections in Section 5. 



3 Results on the leading IW function ^Ai^) 

At leading order, we have extended the SR method to the baryon IW function C,a{'w) 
of the transition A;, — )■ Acii'e [9]. Defining the slope and the curvature of the elastic 
IW function ^a(w) as 

uH = i-pl{w-i) + ^{w-iy + ..., (17) 

we recovered the Bjorken SR for this transition |10j : 



pi 



n>0 
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The quantities t[^\i) denote the = O"*" — t- 1~ IW functions at zero recoil (the 
hght cloud has quantum numbers and n is the radial quantum number). Let us 
point out that only the intermediate states A^"'' with isospin 1 = can contribute 
to the SR. Therefore, Eqn. (ITSl) implies the lower bound for the slope [TO] : 

pi > . (19) 

Using the whole set of SR obtained within the non-forward amplitude method, 
we also obtained a new lower bound on the curvature a\: 

> ^ K + {Plf] ■ (20) 
The bound (l20ll arises from the SR deduced in [9]: 

m = ^J2[^t\i)]' , (21) 

n>0 



2j2\rt\l)]\j2^[''\l)rt^'il) = 0, (22) 



n>0 n>0 n>0 

+4 5^r(")(l)r(")'(l)+a(l) = l. (23) 

n>0 

Eliminating the unknown quantity ^ r^^"'' (l)r^^"^'(l), we finally obtained for the 

n>0 

curvature: 

-I - m = \\pI+ ipir + E [^^'^^'(i)]'} (24) 



5 

that implies the improved bound fl20l) . 

In [15], a group theoretical method to study IW functions has been exposed, that 
sheds a new light on these results. In this approach, a current matrix element splits 
into a heavy quark matrix element and an overlap of the initial and final clouds, 
related to the IW functions, that contain the long distance physics |E]- The light 
cloud belongs to the Hilbert space of a unitary representation of the Lorentz group. 
Decomposing into irreducible representations, one obtains the IW function as an 
integral formula, superposition of irreducible IW functions with positive measures, 
providing positivity bounds on their derivatives. 
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One demonstrated in [15] that this Lorentz group method is equivalent to the SR 
approach, and summarized all the possible constraints. The general formalism was 
thoroughly applied to the case j = for the light cloud, relevant to A5 — Aciue- One 
recovers the bounds f|T9|) and fl20l) and gets new bounds for higher derivatives. Also, 
the Lorentz group approach provides rigorous criteria to decide if a given ansatz 
for the IW function is compatible or not with the general SR, and therefore with 
HQET. 

Let us now recall a result that will become relevant in the present paper. In [15] , 
it has been demonstrated that the k-th derivative of the elastic IW function is given 
by the expectation value of a polynomial of degree k: 

k 



2fc. 



k\ 



n( 



<\\{X + l^ 

i=l 



(2/e + l)! 

with the expectation value defined by 

/•oo 

< f{x) > = f{x) du{x) 



l)> 



(25) 



(26) 



where z/ is a normalized positive measure with support in [0, 00 [ 
One obtains for the first two derivatives: 

1 



pi 



< x> 
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< x{x + 3) > 



(27) 



The point is the following. If the first derivative p\ attains its lowest possible 
value f|T9|l . then one gets for the first moment: 



<x> 



X du{x) 



that completely determines the measure i^: 



< X > = 



du{x) = 6{x) dx 



This implies, for the different moments: 

< x'' > = 



and therefore 



pi =0 



ef(i) = o 



(28) 



(29) 



(30) 



(31) 



As a particular case one gets the result concerning the curvature: 

pi =0 => al =0 (32) 

that we will use below. 

4 SR for the subleading 1/mg form factor A{w) 



From the definition (|T3|) . one gets, inserting intermediate states and taking into 



account that Cmag does not contribute: 



where Ai?*^"^ = m.{n) — m. is the mass difference. Since, as argued above, only 

Ac Ac 

the operator O^'^^ ^ contributes in this relation and Ac"'' is an excited state with the 
same quantum numbers Jj' = |g as the ground state. The only possibility is for 
j = L = 0, which means that in the sum (!33|) only the radial excitations contribute. 

From ( 1331) . we realize that, as it should, ^4(1) = because = 5n,o and the 

slope of A{w) at w; = 1 is given by the expression 

This relation will be used below to obtain a bound on the 0(l/mQ) correction at 
zero recoil of the axial form factor Gi(l). 



5 Sum rule for the 6\ , correction 

Considering the spatial component of the axial-vector current and keeping terms 
of ©(l/mg), one can write a SR similar to the one formulated in the meson case 
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(formulas (114) of [I6] and (5.6) of [IT]): 



G'i(l) 



4ee 



{A^:'\0+,l^){v,s')\A\A,{v,s)) 



s,s' n^O 



Am 



1 



2mr 2my 



3 2mr 2mh 



X 



(35) 



where the parameter —A is the mean kinetic energy value defined by 

1 



A 



2mA, 



■{A,iv)\h^^\tDfh^^^\Mv)) 



(36) 



The factor | in front of the sum (1351) comes from average over the spins s of the 
initial state, and is consistent with the appearance of 1(71(1)1^ in the l.h.s. and also 
with the OPE in the r.h.s.. The r.h.s. of fl35|) is reminiscent of the one found in the 
meson case with the difference that the magnetic part vanishes for the Aq ground 
state Notice that in formula (1351) and in what follows, a single spatial component 
of the axial-vector current should be understood to contribute. 

A similar SR has also been obtained in the baryon case for the vector current 
and with intermediate states = 1~ (J^ ~ I ' I ) P^ - 

In formula (l35l) , one must keep in mind that the matrix element < A^"'' | A | Af, > 
implicitly contains the double insertions of the kinetic and magnetic parts of the 
Lagrangian perturbations to the axial-vector current A'^ = c'-f'^'-f^b. Only the spatial 
part A = C7756 contributes since in the heavy quark limit the component 07^756 
is a subleading corrections to the current itself. Also, in the sum of the l.h.s. of 
(135|) . only the states with the quantum numbers A^J^\j^ = 0+, 1+) contribute, as 
will become clear below. 

Writing Gi{l) under the form: 



Giil)=r]A + S 



(Gi) 
1/ml 



(37) 



we get the following expression for the ^/rriQ corrections: 



-CI 



s,s' n^O 



2mc 2nib 



(Ai"^(0+)(t;,sO I A| Afc(t;,g)) 



8 1 1 
3 2mf. 2mh 



A 



We now consider separately the final states Ac"'*(0"^, l"*") i.e. with quantum num- 
bers = O"*", = i"*" that are attained by Ckin,v insertions, and = 1+, = 
^ and I that are coupled to Cmag,v insertions. 

5.1 Matrix elements (A^^^I^IA^) of transitions 0+ 0+ 

Let us now compute the contribution 

\{A^J'\0'-){v,s')\A\A,{v,s)) p 

4m,(„)m. ^ ^ 

s,s' n^O ' Ai, 

to Eqn. (l38l) . In this relation, the matrix element {A^J^^^ \hi^^ Ah^\o)\Ab) , that con- 
tains the Ckin,v insertions on the initial and final legs, is given by: 

Oil/mQ) 

AE<^[^r^J ^4,n,.„,m,,y,;,„,„;^ u,,MZu.,M. (40) 

The proof of this last relation is as follows. To compute this quantity, we must 
insert the l/mq Lagrangian perturbations on the b and c quark legs. Since we have 
seen that Crnag,v does not contribute when the light cloud has = 0"*", we have (we 
make lighter the formulas by skipping the normalization of the states): 

{A^J'\v,s')\hi'^Ah^^\0)\A,{v,s)) = 

Oil/mg) 

^(A(")(t;,.') \ tJd^xT{Oit^S^),hi'^M\0)} I A,iv,s)) 
1 

2m[, 

Inserting intermediate states, one has: 

{A^:'\v,s') I zJd'xT{Oit^^{x),h^:^Ah^^\0)} I A,{v,s)) = 
-Ak^{^c''\v,s') I Oit^M I A^{v,s')){A^{v,s') I hl^^Ah^^\0) \ A,{v,s)) (42) 

and 

(Ai"Ht;,s) \zJd'xT{O'mx)X''^Ah'^\0)} \ A,{v,s)) = 
+ ^(Ai")(i;,.') I hi^^Ah^^\0) I A^-\v,s)){A^:\v,s) \ O^Sn,M \ A,{v,s)) (43) 
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+ ^(A^J'\v,s') I tJd'xT{O'iS^^^{x)X'^Ah^:'\0)} I A,{v,s)) . (41) 



because Ofj)^ conserves the spin projection. Using heavy quark flavor-spin symmetry, 
the fact that the matrix elements are at zero recoil and that the matrix elements of 
0^2] independent of the spin projection: 

{A,{v, s') I h^:^Ah^:>\0) I A,(^, s)) = (A(")(i;, s') \ h^Hh^\Q) \ A["V, s)) (44) 
and from the relation for the free current (without C insertions): 

{K{v,s')\h^^Ah^\Q)\K,{v,s)) ^ = u\^{v,s')tuK,{v.s), (45) 

O({l/mQ)0) 



the formula (140|) follows. 

5.2 Matrix elements {Kt'^\A\Kb) of transitions 0+ 1+ 

The current A involves a change of flavor b ^ c and the spin matrix element 
u\^{v, s')TjU\^{v, s), as we have seen above. To compute the matrix elements {A^\l~^)\A\Ah), 
we need to take into account that the spin of the final states can be J = | or J = | 
and consider the sums of products of matrix elements (for example, the current A^ 
and the spin projection = +|): 

^(A(")(l+)(i,H-i) I I A,(0+)(i,+i))(A,(0+)(|,+i) I E^T^" I A,(0+)(i,+i)) 



-^^(Ai"ni+)(iH-|) I ^'r^' I Al'^^(l+)(i+|))(Ai")(l+)(i,+i) I OW^,, I A,(0+)(i,+i)) 

(46) 

and 

-^^(Ai"Hl+)(iH-|) I SV'^" I A['^)(l+)(i,+i))(Ai")(l+)(i,+i) I I A,(0+)(i,+i)) 

(47) 

where AE\"''^ = m, („),_,,,— m , Taking into account that the current operator 

^ '(1+) Ac(0+) ^ 

acts on the heavy quark, applying the Wigner-Eckart theorem and summing 
incoherently the contributions with final spin | and |, we obtain an expression 
proportional to 

^ 1 (Ai"^(l+) |1| A,) 1^ 
11 



n 



'a 



(n) 



1 1 



3 2mb 2mc 



2 II ^ II 2 
(n) 



I)' 



1 



1 



1 2V2' 
2mb 3 

1 1 



3 2mc 2mb 



> 



(48) 



5.3 Comparison with the 1/mg corrections to G'i(l) 

Using formulas (4.6), (4.19) and (4.20) from ref. [12], the axial-vector form factor 
inh at zero recoil, Gi(l), can be cast into the form: 



Gi(l) =r/A + i 



X 



+ 



\ 2mc 2mi, J ' 3 2mc 2ini, 
2 y2m. 2m, / ^ 'Dii^) + ^D,{1)] + ^2^4^2(1) (49) 

where the functions (z = 1, 2) correspond to the double insertions of the operators 
O^^n V pTj) and Omag,v (112!) (Q = c, 6) ou the initial and final heavy quark legs (the 
notation Di stands for double insertions). 

Taking the square of ( 149|) and neglecting cross-terms between radiative r]A_ and 
l/mg corrections, one finds, comparing to the SR (l35l) . the relation: 

^[(Ai"H0+,l+)(^) Ml A6(t;)) P _ 

a(") 



4m^(„)m^^ 



1 

2m>; 



[Di(l)- 3^2(1)] -^^8/^2(1). (50) 

Prom poj) and pHl) . we realize that indeed, as expected for coherence, both the l.h.s. 
and r.h.s. of ( l50l) have the same functional dependence in and 



6 Bound on the 0(l/mg) correction at zero recoil 



From the positivity of 



l/ml 



> 



we obtain, from the expression 

2 o 1 1 



the inequality: 



1 



2mr 2mh 



3 2mr 2mh 



A 



s,s' n^O 



{Ai^\0+){v,s')\A\M{v,s)) 



(51) 



12 



that gives, from (HOj) : 



2mr 2my 



1 1 



3 2mc 2mb 



1 / 1 



> -V 

2 y 2mc 2m;, / 



where we have used the equahty 



1 (Ai"n^,«)|C?lI„(0)|Ac(^,^)) 



-1 2 



2 ^ 4mA,mA, 



since it is understood that a single component of S contributes. 
Using now the Schwarz inequahty 



< 



\ 



V n / \ n / 



to Eqn. ([33]): 



A 



W 



we therefore obtain: 



E 



1 (Ai"Ht;,s)|OL^iJ0)|Ae(^,.)) 



1 1 



2mc 2mfe / 3 2mc 2mb 



+ 



2 \2mc 2mb 



[A{w)f 



(52) 



(53) 



(54) 



(55) 



(56) 



Since the inequahty fl56p is vahd for any value of the variable w, we can consider its 
limit for w ^ 1, taking into account that A{1) = and = Sn,Q'- 

2 



-s 



(Gi) 



> 



1 



1 



1 1 



1 / 1 



2mc 2mh ) 3 2mc 2mb 

2 r.//iNi2 



A 



2 \2mc 2mb 



[A'(l)]^ 



'(1) 



2 • 



(57) 
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Prom Eqn. ( 12^ obtained in ref. [9]: 



one gets: 



{n)n 



2 \2 



- Pa - (Pa) 



1 



1 



+ 



2mc 
1 



1 1 



2mb J 3 2mc 2m;, 

2 



A 



1 



10 V2m, 



2mf 



[A'(l)]^ 



(5J 



(59) 



This inequality is one of the main results of the present paper. 

From (IT^ and 0161) . we see that the l/mg correction at zero recoil —^^^2 is 
crucial in the extrapolation of the differential rate (ir(A;, — )■ h.i.tP^)/dw near zero 
recoil. This is needed to check that the value of \Vcb\ that would fit the exclusive 
baryon semileptonic data is indeed consistent with what we presently know on this 
parameter from the meson exclusive and inclusive determinations. It is in this 
respect that the bound (!59|l has to be used. 

To get a quantitative estimate on the bound (1591) . we need information on the 
parameters —A, p\, a\ and A'{1). These quantities are on quite different footing 
and their plausible values will be discussed in Sections 9, 10 and 11. 

7 Correlation between v4'(l) and the shape of the 
leading IW function ^k{w) 



Since the inequality fl59|) holds for any values of p\ and a\ satisfying the general 
constraints ( IT9l) and ( l20l) . it should hold, in particular, for their lowest values. In- 
deed, one has that if p\ and a\ attain their lowest values, the r.h.s. of the precedent 
inequality would diverge. This feature seems completely unphysical. Therefore, we 
expect a correlation between A'(l) and a^, namely, one should expect that A'{1) 
vanishes if (j\ attains its lowest bound ( 1201) : 

A'(1)^0 for al^l[pl + {plf]- (60) 

14 



We give now an explicit and independent demonstration of this interesting fea- 



ture. 



Let us consider the derivative of ( 133|) at zero recoil : 

1 ,(„y„,(Ai"'(».s)|Olt„(0)|M».s)> 



Using again the Schwarz inequality as in Section 6, we obtain: 

'(A^")(i;,.)|Oit,(0)|A,(t.,s))' 



and from relation fl24l) we obtain: 



4m,{n)mA,,/f°(„a'° 



[A\i)r<lU-l[pi+{pi)]W 



Therefore, if the lower bound fl20|) is saturated, one obtains 



=0 



(61) 



(62) 



(63) 



(64) 



as we see from inspection of relation flB5]) . 

On the other hand, one has the important result (132|) . i.e. p\ =0 implies 
a\ = 0. This feature has a clear theoretical-group interpretation [15], that implies: 

pi = =J> v4'(l) = . (65) 

Formulas fl64p and (1651) are quite nontrivial because they imply a strong corre- 
lation between the shapes of the leading IW function iA{ui) and of the subleading 
form factor A{w). 



8 Parameters involved in the inequality on —Sy^2 

Q 

The inequality (!59l) involves, besides the quark masses rric and nib, the slope p\ and 
the curvature of the elastic leading IW function C,a{w), the mean kinetic energy 
value —A and A'{1), the derivative of the subleading form factor A{w) (IT5]) at 
zero recoil. In order to compute the r.h.s. of fl59l) . we need an estimate of these 
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quantities. The theoretical work on the IW function ^a(w) and the subleading form 
factor A{w) have not been so exhaustive as the corresponding ones in B meson 
semileptonic decays. We will now make a review of what can be known about these 
parameters, namely, within HQET, the QCDSR approach and the quark model. 
Notice that we concentrate on the baryon Aq and leave aside other heavy baryons. 



9 Discussion of the parameter —A 



Let us first discuss the parameter —A. In principle, one should use its measured 
value extracted from a fit to the data on the semileptonic inclusive decay A;, — )■ Xciue, 
as has been done in the case of mesons. However, to our knowledge, this analysis is 
not available yet. For the time being, following |1T], we will use a lower bound on 
—A that can be obtained from data on the heavy baryon spectrum. 

Let us write down the SR ( 1T8|) for the slope p\ and the first and second moments 
of this SR [H], that give the following expressions for A and —A: 

2 



pI 



n>0 

A = 2 ^(A'^'^) - A) [rf^(l) 

n>0 

-A = 3 ^(A'(") - Kf 

n.>0 '" 

At leading order in 1 / niQ , A and A'*^"^ are given by 



(66) 
(67) 
(68) 



A 

A'(") 



"^Aq(0+) - , 



(69) 
(70) 



As pointed out in [TT], Eqns. (l67|l and fl68|) can be combined to give the bound: 



A > -A(A' - A) 



(71) 



where A' = A'*^''^ i.e. A'*^"^ for n = 0. The inequality (17T!) . that will be used to 
bound the OPE term in the r.h.s. of the inequality 0591) . follows from the fact that 
A'(") — A > A'*^*^) — A for 72 7^ 0, the n = state assumed to be the lowest mass state, 
a very reasonable hypothesis. 
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In ref. [TT], A has been estimated from the masses of the = O"*", = \ 
states Ac(2.286), Afe(5.620), and A' — A from the hehcity weighted average for the 
= 1" doublet Ac(|~)(2.595), Ac(|~)(2.625). The following results were obtained 

A ~ 0.8 GeV , A' - A ~ 0.2 GeV . (72) 

Then, from the inequality ( ffTl) . one obtains : 

- A > 0.24 GeV^ . (73) 



10 Theoretical results on ^a(^) and A(w) 
10.1 Heavy Quark Effective Theory 

We have already commented the basic work done within HQET for the exclusive 
decay A;, — )■ A^iVi. Using the HQET formulation of the heavy hadron fields of 
arbitrary spin |6] , one has recently extended non-forward SR method to the baryon 
IW function ^a(w) [9], recovering the lower bound for the slope p\ = —^^(1) > 
[To], and generalizing it for higher derivatives. In particular, the new improved 



bound ( 120|) on the curvature > | [p\ + {PaY] has been derived. 

With the method based on the Lorentz group [15], one has formulated a general 
criterium to verify if a given phenomenological ansatz for the IW function C,a{w) 
is consistent with the sum rules, that in turn are equivalent to the Lorentz group 
approach. 

Among the various phenomenological results listed in [15], it will be useful for 
our purpose here to recall two possible parametrizations for ^a{w). On the one hand, 
the "dipole" form: 



,2 



aH=(^r' (74) 



is fully consistent, provided the slope satisfies 

1 



The curvature is then given by 



pi > 4 • (75) 



^1 = f + ipir ■ (76) 
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On the other hand, the form {w = ch(r)): 



sh ^Ta/I — 3p 



is consistent for the whole range of the slope 

Pa > (78) 



allowing saturation of the general lower bound for the slope obtained in [10] . In this 
latter case, the curvature is given by 

i.e. the lower bound fl20l) is saturated and, according to fl64l) . A'{1) vanishes. Of 
course, this form is an extreme case, constructed in [15] in order to have saturation 
of the lower bound. 

We will later use the form (!74l) for ^tv{w) in order to estimated the lower bound 
on —6: 1 2 • 

As far as the subleading form factor A{w) is concerned, HQET implies Eqn. 
(|T^ ^(1) = 0, and the present paper has established two new results, namely, 
^\ ~ ffpi + (Pa)^] ^ ^'(1) = (IMl) . and with the Lorentz group approach, 
p2 = ^ A'(l) = (EHl). 



10.2 QCD Sum Rules 

As it is well-known, unlike the case of mesons, there are two possible interpolating 
baryon currents j\ and j^. 

The simplest hadronic quantities that have been computed within the QCDSR 
approach are the Ag baryon bound state energy A and the decay constant /a, that, 
at leading order, are independent of the choice of the current. 

10.2.1 The value of A 

Concerning A for the baryon Ag, there is consistency between the different calcu- 
lations. At leading order, one obtains A = 0.78 GeV in the work of Grozin and 
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Yakovlev [H], and A = (0.79 ± 0.05) GeV in the one of Dai et al. [I9], later con- 
firmed in ref. [20]. In a series of papers by Groote et al. |21], tlie Af, baryon mass 
and decay constant with the radiative corrections at NLO in have been studied, 
with the finding A ~ 0.780 GeV. 

10.2.2 The decay constant /a 

Concerning the decay constant /a, since the pioneering work of Shuryak [22], the 
value obtained has somewhat evolved over the years. In ref. [18], the range /a = 
(1.8 — 2.7) X 10"^ GeV^ was proposed. In the works of Dai et al. [I9] [23], the value for 
its square was obtained /I = (0.3±0.1) x 10^^ GeV^ {i.e. /a ~ (1.7±0.3) x 10^ GeV), 
confirmed in a much later paper [20]. In the more recent work [2l], a somewhat 
larger value is found, roughly = (1.3 ± 0.7) x 10^^ GeV^. In the papers by 
Groote et al. [21], taking into account QCD radiative corrections at NLO, one finds 
/a ~ 0.028 GeV^ that is /| ^ 0.78 x 10-^ GeV^ i.e. a larger value by roughly a 
factor of two than that at leading order. 

10.2.3 The leading Isgur-Wise function ^h{w) 

Let us shift now to the three-point function, i.e. the determination of the IW 
function ^\{w). Dai et al. [23] performed a QCDSR analysis using both currents H 
and ^2 and they found for the slope p\ = 0.55 ± 0.15. Huang et al. [21] reanalyzed 
within the QCDSR the function C,a{'w) and found a higher slope, p\ = 1.35 ± 0.12. 

10.2.4 Power corrections and A'{1) from QCDSR 

To get an estimate of the r.h.s. of the inequality fl59l) . we need the quantity A'{1). 
Let us now discuss the literature on the subleading form factor A{w). 

In HQET, as we have seen above, the seminal paper is the analysis of the power 
corrections in HQET of baryon form factors [12]. The kinetic part of the Lagrangian 
perturbation at 0{l/mQ) is parametrized by A{w) according to ( fT3l) . that satisfies 
the condition at zero recoil A{1) = ( fT4l) . In the present paper, we have seen that 
there are other rigorous conditions on A'{1) from HQET, namely flMl) and fl65l) . 

We turn to the QCDSR estimate of the function A{w). To our knowledge, 
there is only the paper by Dai et al. [23] that has computed A{w), given in terms 
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of a single function J{w) besides the leading IW function ^Aiw), with A{w) oc 
J{w) — C,a{w)J{1), satisfying then A{1) = 0. 

We find here the following question. We have demonstrated above that within 
HQET, the slope A'{1) must satisfy the conditions flMl) and fl65|) that relate the 
leading and subleading IW functions, at least formally. Formally, because of course 
in the real QCD, p\ needs not to be zero. 

Within the QCDSR approach, according to the generic relation, A'{1) depends 
on the slope p\, but it is independent of the curvature a\. Therefore, from the two 
conditions ( IMl) and (l65l) . we can ask whether the QCDSR calculation satisfies, at 
least in a formal way, the second condition ( 1651) . p\ = ^ A'{1) = 0. 

As far as the case of the current is concerned, for p\ = all the terms in the 
corresponding expression for A'{1) are positive-definite and therefore the condition 
(l65l) cannot be fulfilled. On the contrary, for the case of the current jg, for p\ = 
the vanishing of A'{1) is possible, although not in the stability plateau. 

By the way, notice that the condition p\ = ^ A'{1) = could provide a 
consistency check for choosing the right baryonic current, that happens to be J2. 
Interestingly, this constraint seems to lift the ambiguity in the choice of the current. 

However, choosing the current we have obtained lower bounds, depending on 
the IW function slope, that are large, leading to a lower bound on — (5^*f^2 that 
seems to us much too large, as it can for some cases approach 50 %. 

10.2.5 Summary on QCDSR 

For a number of reasons, we do not proceed any further with the QCDSR method 
to extract A'{1). First, although the method gives very stable values for A, the 
values obtained for the decay constant /a do not show such great stability. Notice 
that this constant appears also in the determination of the three-point function. 
Furthermore, as we have seen, f\ depends strongly on the radiative corrections. On 
the other hand, we have seen that the study and results of the elastic leading IW 
function $,a{'w) have evolved with time. Concerning A'{1), the quantity of interest 
to us as far as the bound ( l59l) is concerned, there is the paper [23] that allows 
to compute it. However, there are subtle points that remain to be investigated, 
in particular, the fulfillement by the QCDSR approach of the HQET condition. 



20 



p2 = o^A'(i) = §B- 

In view of the involvement of the study of the QCDSR method for the purpose of 
the present paper, we leave it for future work, and turn now to a simpler approach, 
namely the quark model. As we will see, it allows to compute a value for A'[l) that 
furthermore satisfies the condition ( !65!) . 

11 Calculation of A'[l) in the quark model 

11.1 General formulas 

Let us consider the three four-vectors: 

V = {v°,0,0,v^) , v' = {v'°,0,0,v") and n= (0,1,0,0) (80) 

satisfying v"^ = v'"^ = 1, = —1 and v ■ n = v' ■ n = 0. Then, from the matrix 
element ([1]), we obtain: 

{A,{v', s') I c^b I Ab{v, s)) = (Fi + F2 + Fsw) ua^{v', s')uaSv, s) , (81) 
(A,(t;', s') I cfb I Ab{v, s)) = (Fi + F2W + F3) UA,(t/, s')uA,{v, s) , (82) 
{Ac{v', s') I cfib I Ab{v, s)) = Fi ua,{v', s)iiUA^{v, s) (83) 

where we have used the constraints '/)U\^{v, s) = ua^{v, s) and ua^{v', s')f = ua^{v', s') 
It is very convenient to perform the calculation in a particular frame, namely 
the equal velocity frame (EVF), since in this frame the calculation is symmetric in 
the exchange 6 -H- c. In this frame, one has: 

q = PAc- Pa, , Pa, = ^ q , Pa, = . q 84 

implying = —v'^, = v'^ and 

q' = {rriA. + mA,r-^ . (85) 

In what follows, we take spinors for the bound states as well as for the quarks 
normalized in the following way: 



u{p, s) = \l — — I I with uu = 1 . (86) 



2m \ ^-p 
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Of course, the final results are independent of the adopted normalization. We find 
in the EVF some formulas that will be useful below: 



E, 



At, 



^^^^ and ua^{v',s')ua^{v,s) = {xl^XA,) ■ 



(87) 

Prom the preceding formulas, we find after some spinor algebra at the bound state 
level, using ( 185|) and particularizing to suitable spin projections: 



Fi{w) 
F2{w) 



w - 1 



{A,{v',+l)\cfib\A,iv,-'^)) 



1)1 1 



w + 1 

+ 



w + 1 w — 1 



{A,iv',+'^)\c{-^ + wf)b\Ab{v,+l)) 



w — 1 



{A,{v\+l)\cfib\Ak{v,-l)) 



^9) 



F.(w) 



1 



w + 1 
+ 



w + 1 w — 1 



{A,{v',+l)\c{w^-f)b\At,{v,+l)) 



w — 1 



(A,(t;',+i) |c/i6| A6(t;,-i)) 



(90) 



As we will check below, after computing within the quark model the matrix 
elements present in the r.h.s. of fl88l) - fl90l) . the singularities in and in -^7=^ will 
disappear. 



11.2 Calculation in the quark model 

We now perform the calculation in the quark model of the form factors Fi{w) {i = 
1,2,3) ([3]), (E]) and (E]), together with ((Tj) and (IE]). As we will see, the results fulfill 
this form expected from HQET, and this will allow us to extract A'{1). 

Let us consider the harmonic-oscillator Hamiltonian for baryons {i = 1,2,3): 



I Kj 



(91) 



where K is independent of flavor. Notice that we use this simple potential because it 
is the only confining one that allows to separate the center-of-mass in the three-body 
case of baryons. 
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From now on, we consider the system in which we are here interested, where 
mi = m2 = m (hght quark masses) and ms = mq (heavy quark mass, Q = b,c). 
For details on what follows on the quark model wave functions, see Appendix A of 
1251. The total and relative coordinates read 



R 



p = -^(ri - raj 



and the corresponding conjugate momenta are 

1 



PR = ^Pi 



Pp = -^i.Pi-P2) 



Px 



A = -^(fi + f2-2r3) , (92) 



rriQ ^ 2m ^ 

{P1+P2) - ^ P3 



In terms of these variables the Hamiltonian reads 



H 



Pr 



Pp 



+ 



Px 



\2m+mQ I 



(93) 



(94) 



In the calculations that follow, we will use the wave functions in momentum 
space. The Schrodinger equation gives, for the ground state, the internal wave 
function: 



with 



'^{Pp.Px) = Noexp 



4 2m + niQ 



'^ip'pR'p+plRl) 



(95) 



3 d3 \ 2 



3mr 



Nn 



f 3V3Rf,R^ 



(96) 



We will now calculate in the quark model the matrix elements in (I88l)-( l90l) : 



(Ae(i;',+i) |c)^6| A,(^;,+i)) 
{A,{v',+l)\cfb\A,iv,+l)) 
(A,(^/,+i)|c/iHA,(^^-|)) 



{A,{v', +1) I c (7V - b I Ab{v, +1)) , (97) 
(A,(^;', +1) I c (7V0 - 7V^) b I A,{v, +|)) , (98) 
-(A,(t/,+|)|c7^'6|A,(t;,-l)). (99) 



We need therefore to compute the matrix elements (Ac(w', +|) | c-f^b \ Af,(i;, +|)), 
(Ac(w',+|) I cYb I Ah{v,+l)) and (Ac(w', +|) | c^b \ Ab{v,-^)). Using the nota- 
tion r = {7^,7^,7^}, we have to compute in the quark model, denoting by 3 the 
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mA,, 



active heavy quark: 

{ip^"'^ I u^''\p^)Tu^''\p3) I ?/'(^)) = J dpidp2dp3dp[dp2dp'3 5{pi - p[) 5(^2-^2) 

xS{p[+P^+P3 - ^r.+m^, ^) ^{P^ +P2+P3 
X {pip^,P3) n(^)(p;)rnW(p3) ^("^ (P1,P2,P3) . 

(100) 

We obtain, after some algebra and change of variables: 

(^(^) I u(^\P3)Tu(''\p,) I ^(^)) = dp,dp,^^^^^{p,,p, - 

(101) 

Some important words of caution are in order here concerning the terms that we 
keep in the non-relativistic expansion. Following the usual methods in the quark 
model, we expand the quark spinor matrix element in (11 01 1) 

^^'^ (-\fipx + ^^gVnW f-./|pA - ^^q] (102) 

up to first order in l/mq and then we perform the integral on the internal variables. 

On the other hand, we keep terms of the order i?^m^, that are of order (|) ^, 
corresponding to the non-relativistic limit of the slope of the IW function (for a 
discussion, see [26]). An important remark to be made also is that the HQET 
parameter A can be decomposed, in the quark model, under the form A = 2m + e, 
where m is the constituent light quark mass and e is the binding energy in the 
non-relativistic quark model. Of course, both cannot be distinguished in HQET. 
However, in a non-relativistic expansion as adopted here, e is of order -^5^ that, 
relatively to the constituent mass m is a subleading quantity ^ ~ (f)^- Therefore, 
from now on, we will neglect e and take 2m = A. 

On the other hand, to extract A'{1) we need to go to the orders w — 1, ^ and 
^(tf — 1), that we will also keep in our calculation. 

To summarize, in our quark model calculation we keep the following orders in 
the non-relativistic expansion: leading order in HQET ( ^ j (keeping the order 
A^R^ and the first order in w — 1), order — and order —(w — 1). 
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A lengthy and careful calculation gives the following results: 



(A,(i;',+i)|c7%|Ab(t;,+i)) = FH, 



(103) 



w -1 / A A 



2 V 2mfe 2mc 



w - 1 



(A,(t/, +1) I C7^'6 I A,{v, -i)) = J^^^ f 1 ^ 



where 



F{w) = exp 



^ A^i?^ I 1 + 



A A 



4 \ Ami, Arric 



iw-1] 



(104) 
(105) 

(106) 



i.e. we obtain in the heavy quark limit the IW function F[w) — > C,{w) and its slope 

^2. 



^a{w) = exp 



4 ^ 



with pI = ^ A'R' . (107) 



One has demonstrated elsewhere [15] that, strictly speaking, the exponential form 
( I107p is not consistent with Lorentz group constraints. However, not unexpectedly 
in the non-relativistic harmonic-oscillator quark model, one has such a behaviour 
in the EVF. However, this is not really an inconvenience since the non-relativistic 
expansion cannot hold at large w — 1, but only near zero recoil, in the neighborhood 
of w = 1. Therefore, we will only consider the first values in the Taylor expansion 
of ^a{w) in f llOTp . namely the normalization ^a(I) = 1 and the slope —^^(1) = p\. 
Interestingly, the value i?^ ^ 6 GeV~^ extracted from the light baryon spectrum [25] 
gives numerically a value for the slope that is not non-sense (A ~ 2m ~ 0.70 GeV): 



pI ^ 1.27 . 



(108) 



We now go back to the calculation. Taking into account that, from (18^ and 
51): 



V = —V 



w — 1 



and 



= 



w + 1 



(109) 
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we can now compute (1971) - (!99l) : 



{A,{v', + '^)\c^b\Ai,{v,+l)) = F{w) 



{A,{v',+l)\cfb\A,{v,+'^)) = F{w) 



w 



1 w-1 / A A 



2mb 2mc 



w + 1 w -1 / A A 
+ —r- 



2mh 2m, 



{A,{v\+l) I cfib I A,(t;,-|)) = -F{w)^ 



w -I ( A A 
1 + - — + 



2mh 2m, 



w — 1 
4 

(110) 

w — 1 
4 

(111) 
(112) 



Inserting these expressions in the formulas (155]) - (M?]) . we obtain: 



F^{w) = F{w) ( 1 + + ^ 



2mh 2mc 



(113) 



Foiw) 



F.(w) 



F{w) 
w + 1 



F{w) 
w + 1 



w + 1 / A A 



2mh 2m, 



w — 1 



A A 



w + 1 f A A 



2m,. 2mh 



w — \ 



2mh 2m, 



A A 



2mh 2m, 



;ii4) 



;ii5) 



It is important to point out that the apparent singularities in the expressions (158)) - 
(190|) have disappeared, since the singularities at the bound state level are cancelled 
by their inverses at the quark level. Notice the different role of mb and in the 
first term of the expressions for F2{w) ( 1114!) and F^{w) ( 11151) . 



11.3 Final results and comparison with HQET 

From the expressions (I113p - (lll5p . we can compute the quantities Fi{l) {i = 1, 2, 3): 

A A 



^2(1) 



2mh 2mr 



A 



2mc ' 

A 
2mb 



116) 
117) 
118) 
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and the derivatives (i = 1,2,3): 

^2(1) 

1 A „ 

^3(1) = oT^ + P' 



-p 

1 A 
22m^ 
1 A 









H 


( ^ H 






\2mh 


2mJ _ 



p 



A 

2mc 

A 



(119) 
(120) 
(121) 



2 2mb 2mf, 

Let us now compare with the HQET formulation, Eqns. ([3]), ([5]) and (Q. We 
first realize that there is consistency between the general structure in ^ of F2{w) 
and Fslw) and the quark model results. 

In the quark model, the function B2{w) (I8j) satisfies: 

A 



B2{1) 



-A and 



i?2(l) = ^ + /A- 



;i22) 



Notice the very non-trivial point that -82(1); the derivative B2{1) (11221) as well as 
Fi{l) (11161) coincide exactly with the HQET results, as can be read from the general 
formulas ([3]), (ED and ([8]). 

Let us now compare the first derivative F[{1) between the quark model and 
HQET. In HQET, one has from (JHl): 

1 1 



-P 



2mh 2mc 
that can be expressed using ([7]) in the form 



Fiil) 



-P 



1 



2mh 2mr 



[B[{1) - B'^il)] 



^ + A'(l) - 5^(1) 



(123) 



(124) 



with -62(1) given by (I122p . We now compare with the quark model result (I119p . 
This yields the equality: 
1 1 



P 



2mh 2771. 



A ,,, , A 







f A 




= V 


H 
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2mc/ 



from which can be extracted the expression for A'{1): 

A . 



A'il] 



-P 



(125) 



(126) 



which satisfies the constraint A'{1) — )■ for — )■ and which is the main result of 
this Section. 
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12 Lower bound on — 

It might be enlightening to spht the inequahty flM?]) into two individual contributions, 
one from OPE, the other from the sum of intermediate states (IS), respectively: 

- > AoPE + A,5 (127) 



where 



1 

^OPE --- 



2m, 2m/, / 3 2mc 2m;, 



A (128) 



and 

/V _ 3 / 1 1 y [A'(i)]2 

- To ^2;;^ " 2^;^; ai-|[pi + (pi)2] • ^^^^^ 

With the masses m, = 1.1 GeV, m^ = 4.5 GeV and the bound -A > 0.24 GeV^ 
(I73|) . one gets: 

AoPi? > 0.03 . (130) 

To bound Ajs, one needs not only information on A'{1) but also the shape of 
the IW function to have the slope p\ and the curvature a\. 

For the leading IW function, we adopt here the "dipole" shape (17il) . that implies 
(1761) . For the quantity A'{1), we adopt the expression 01261) obtained in the quark 
model. We use the reasonable value A = 0.7 GeV, twice the light quark constituent 
mass, that agrees qualitatively with the QCDSR calculations, as we have seen in 
Section 10.2.1. Then, A/5 depends only on the slope p^, that we vary. We find: 

(131) 
(132) 
(133) 

Of course, one must keep in mind that the contribution of the l"*" states in the 
sum (138 p has not been estimated, only their positivity is taken into account. This 
suggests that the actual value of — <^J/,^2 could be relatively large. Another remark 
to be made is that in the baryon case, there is no strong cancellation like for the 



for 


pI = 


0.5 , 




> 0.052 


for 


pI = 


1.0 , 




> 0.045 


for 


pI = 


1.5 , 




> 0.043 
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bound that appeared in the meson case due to the various subleading form factors 
Xii'w) (i = 1,2,3) [1]. A single form factor A{w) contributes in the present case. 

To conclude, one gets a relatively large lower bound for —5y^2 , although, of 
course, the final numerical results rely on two phenomenological hypotheses: the 
"dipole" shape ( j74|) for ^{w) and the quark model estimate ( 1126^ for A'{1). 

It is quite important to have an estimation of A'{1) in the QCDSR approach. As 
stated above, one needs to clarify in detail some important points of this theoretical 
method, and we plan to do it in near future. 

The crucial decisive next step would be to have a measurement on the lattice of 
the correction at zero recoil —Sy^2 , as has been done already for the meson case. 

13 Conclusions 

We have obtained new results in the heavy quark expansion of Heavy Quark Effective 
Theory that are relevant for the differential rate of the decay Af, — )■ Acii^i, that will be 
measured in greater precision at LHCb. We have written down a SR for the elastic 
subleading form factor A{w) at order l/rriQ, that originates from the Lagrangian 
perturbation Ckm- This result, together with another SR in the forward direction 
for the axial-vector form factor has allowed us to obtain a lower bound 

for the correction at zero recoil — ^i'f^l in terms of the derivative A'lfl) and the 

l/mg ' 

slope p\ and curvature cr^ of the leading Isgur-Wise function C,a{'w). We have found 
that the derivative A'{1) must vanish in the formal limits a\ ^[{PaY + Pa] and 
p\ — )■ 0, establishing a non-trivial correlation between the shape of the leading IW 
function ^a{w) and the subleading one A{w). We have discussed results of the theory 
for the leading and subleading IW functions (HQET and QCDSR). Moreover, we 
have performed a calculation in the quark model of the leading and subleading form 
factors, that agree with HQET up to order and including [w — 1). Consistently, the 
result for v4'^(l) in the quark model vanishes in the limit p\ ^ 0. We finally bound 
from below the l/m^ correction —6[^^\ to the axial form factor at w = 1, that 

/ Q l/m^ ' 

contributes to the differential rate of A;, — )■ Acit'e at zero recoil. 
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